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Abstract 

This article presents a technical innovation consisting of replacing the piston-cylinder sys¬ 
tem by metallic bellows in free-displacer Stirling engines. It deals with the advantage of heat 
engines of low power, integrating this innovation with plane solar collectors. A theoretical 
study requiring notions of finite time thermodynamics determines the optimised dimension¬ 
ing of such systems and makes them appear serious competitors, on the one hand to solar 
cells to exploit solar energy and, on the other hand, to power stations regarding kW h unit 
cost. 
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1. Introduction 

Energy is the essential factor for higher standards of living. It is, therefore, 
absolutely vital in economically developed societies. Considered in the long term, 
the question of maintaining higher standards of living is intrinsically linked to 
hypothetical sources that can deliver energy in profusion, without a limit in time. It 
is easy to see that solar electromagnetic radiation provides such a source. Using an 
elementary calculation, we can see that the surface of a square of side 400 km, 
equipped to transform solar electromagnetic radiation into electrical energy with 
0.15 efficiency, is sufficient to satisfy the energy requirements of the world, at the 
current level of approximately 10 TW. Thus considered, the question of long-term 
economic development meets the recommendations of the Harare declaration [1], 
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Nomenclature 

A surface area [nr] 

c p specific heat at constant pressure [J/kg K] 

c v specific heat at constant volume [J/kg K] 

h heat exchange coefficient [W/nr K] 

m mass of driving medium [kg] 

p pressure [Pa] 

Q heat exchanged during a cycle [J] 

Q heat power [W] 

r mass constant of gases [J/kg K] 

S entropy [J/K] 

T temperature [K] 

U internal energy [J] 

V volume [m ] 

W work exchanged during a cycle [J] 

W mechanical power [W] 

Greek symbols 

a regenerator effectiveness 

0 flux heat [W/m 2 ] 

r} efficiency 

Indexes 

abs absorber 

amb ambient atmosphere 

h hot source for driving medium 

coll collector 

hot hot source for engine 

c cold source for driving medium 

cold cold source for engine 

lk heat leakages 

fric viscous or static friction 

g global 

M absorber boundary value, see Eq. (9) 
max driving medium in isothermal expansion 
min driving medium in isothermal compression 
mot motor 

opt optimal 

reg regenerator 

S solar 

1, 2, 3, 4 points of the cycle 
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Intermediary variables 

P, Po Eq. (12) 

ao Eq. (19) 

bo, b u b 2 Eq. (21) 

K , c 0 , ci, c 2 Eq. (22) 

a 1? Pi Eq. (25) 

a 2 , P 2 Eq. (26) 


advocating the research, development and implementation of technical systems for 
the massive exploitation of solar energy. 

The value mentioned above of 0.15 corresponds to the average efficiency of solar 
cells, systems generally chosen to exploit solar electromagnetic radiation. The 
efficiency of a cell is, more accurately, 0.12 regarding commercial quality, 0.17 for 
improved quality and 0.24 for the best cells manufactured in laboratories. These 
values, all notably inferior to those of heat engine efficiency used in conventional 
energy production, would nevertheless be satisfactory within the context of exploit¬ 
ing solar energy—exploiting any other form of energy originating from solar 
energy (wind, hydraulic or biomass) does indeed give derisory final efficiency, 
inferior to 0.001—but ultimately, it is essentially for economically competitive rea¬ 
sons with regard to conventional energy production that solar cells are not more 
widely used. Following this line of thinking, the heat option for exploiting solar 
electromagnetic radiation should be taken into consideration as the technology of 
heat engines allows for moderate cost prices and their conception, today, could 
benefit from solid theoretical knowledge (consult the works by Walker [2], Organ 
[3], Senft [4] and Bejan [5]), when combined with sustained creativity (see, for 
example, the patents of Lamos [6], Bliesner [7] and Kinnersly [8]). 

The author of this article would like to make credible the idea of exploiting the 
field of solar electromagnetic radiation using the heat option, by publishing a tech¬ 
nical innovation, apparently entirely original in bibliographies and patent registra¬ 
tions. This innovation consists of replacing the traditional piston-cylinder system in 
the Stirling cycle heat engine with metallic bellows; the main obstacle to achieving 
very high efficiency with heat engines from 100 to 3000 W of power therefore being 
eliminated, it becomes necessary to study what level of optimal performance the 
combination, composed of such a heat engine and a solar collector, could reason¬ 
ably attain. 

The working principle of the system is justified and described in Section 2. The 
mathematical model is set out in Section 3 and optimal dimensioning is determined 
in Section 4. The numerical application in Section 5 enables the most sensitive 
technical points to be defined and tackles, finally, the question of economic compe¬ 
titiveness in terms of the cost of kW h produced. 
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2. Solar heat engines studied 

2.1. Preliminary on Carnot efficiency cycles 

The Carnot cycle is composed of isothermal compression (1-2) at temperature 
Tmin, adiabatic compression (2-3), isothermal expansion (3^1) at temperature T max 
and adiabatic expansion (4-1) (Fig. 1). When it occurs under infinitesimal tempera¬ 
ture differences, the temperatures T max and T min become identified with the tem¬ 
peratures T h and T c of the heat sources, and in the absence of friction, the 
efficiency of the cycle attains the well-known upper boundary value rjcamou depen¬ 
dent only on the temperatures of the hot and cold sources: 


^7 Carnot 

By explaining the terms involved in defining the efficiency of an engine 
rjmot = —W/Qh 0 t, (1) can be established. Using first principles, the reverse of the 
indicated mechanical work, exchanged during a cycle, can be written: 

~W= (2hot + 2cold (2) 




Fig. 1. The Carnot cycle and Stirling cycle. 
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and we can express the efficiency, written according to the heat exchanges with the 
heat sources, as: 


h mot 


Qh.oi T Scold 


Sh 


ot 



Relation (3) then leads Eq. ( 1 ), using second principles: the heat values Shot and 
Scold are replaced, exchanged by the driving medium during a reversible cycle, by 
7h • (S 4 — S 3 ) and r c • (S 2 — Si), respectively, and it can be noted that — (S 2 — Si) = 
S 4 — S 3 in the Carnot cycle. 

The existence of an infinite number of cycles reaching Carnot cycle efficiency, or 
Stirling cycles, is well-known (see, for example, [9]). It can be observed by remark¬ 
ing that all relations of the form 


h mot 


Shot Sreg Scold T Sreg 


Sh 


ot 



replacing Eq. (3) lead to Eq. (1) using the calculation indicated above. The terms 
— Sreg an d Sreg may be considered as the heat exchanged during the evolutions 
(4-1) and (2-3), and verifying by definition, for example: 


3 



TdS = 

2 


*1 

TdS. 


4 



As the Sreg and —Sreg exchanges take place without disturbing the heat exchan¬ 
ges with the hot and cold sources, the heat value Sreg must be stored and des¬ 
tocked in an intermittent manner within an annexed device, the regenerator. 

Finally, Stirling cycles are composed of a compression (1-2) and expansion (3-4), 
both isothermals and verifying — (S 2 — S\) = S 4 — £ 3 , as well as the two inter¬ 
mediary phases (2-3) and (1-4), whose evolution is only limited in that the heat 
| Sreg 1 5 exchanged with the regenerator, verifies (5) (Fig. 1 ). 


2.2. Description of the innovative solar heat engine 


The objective of this article is to determine the optimal performances that a 
competitive solar heat engine can realistically reach, on the energy and economic 
planes. 

Solar electromagnetic radiation is collected using a flat plate solar collector, 
which is easier to produce than a concentrating collector. The heat power, originat¬ 
ing from the solar electromagnetic radiation and transmitted to the driving 
medium through an absorber, is transformed into mechanical power by a heat 
engine, whose theoretical cycle attains Carnot efficiency. 

The succession of isentrops and isotherms in the Carnot cycle requires the con¬ 
ditions of heat exchange to be modified several times per cycle. Technically, this 
leads to complicated systems and, to ensure the reliability and economic competi¬ 
tiveness of the solar heat engine, in this case, a free-displacer Stirling engine was 
chosen. The configuration of a free piston and regenerator was adopted, which, 
according Rochelle [10], is the most appropriate. 
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Fig. 2. Diagram of the innovative solar heat engine. 


The practical arrangements which bring the true efficiency of the heat engine 
close to theoretical efficiency are shown in Fig. 2. The three main points are as fol¬ 
lows: 

1. The piston-cylinder system is replaced by metallic bellows, reducing the mechan¬ 
ical friction to a negligible value, 

2. The regenerator, interdependent with the working piston, has a weak flow area 
with regard to heat exchange surface; it is formed from a stack of air-permeable 
plates, with heat insulation between them, so the storage and destocking of heat 
takes place under the temperature gradient imposed by the absorber and fan 
cooler, 

3. The bases of the cylinder, absorber and fan cooler have a large surface area, so 
heat exchange can take place with only slight differences in temperature. 


3. Mathematical model 

3.1. Hypotheses 

The mathematical model is a synthesis of those by Magot-Cuvru [9] and Badescu 
et al. [11]. Similar to the latter, the evolutions (1-2) and (3-4), taken to be revers¬ 
ible isothermals, take into account the difference in temperature between the sour¬ 
ces and the driving medium. The absorber acts as the hot source; thus, = T^ s . 
The cold source for the engine is the ambient atmosphere; T c = T am b is obtained. 
Finally, in normal working conditions, we obtain T m [ n > r am b and T max < T^ s 
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Ambient Atmosphere 




-cold 


T C = T. 


amb 


Ambient Atmosphere 

Fig. 3. Modelling of energy exchanges. 



(Fig. 1). The driving medium behaves as a perfect gas of sensitive specific heat 
values c v and c p and mass constant of gas r. 

The energy exchanges take place at a stationary rate of flow. The fluxes taken 
into account are indicated in diagrammatic form in Fig. 3. All the fluxes exchanged 
outside the engine are modelled in the same way as by Badescu et al. [11]. The heat 
power Q$, originating from the solar electromagnetic radiation, is redistributed by 

the absorber, for 2hot ? to the engine, and for 2ik ? into the ambient atmosphere, in 
the form of convective heat leakages; we obtain 


2s + 2hot + 2lk — 0- 


(6) 





















1832 


L.-C. Valdes / Renewable Energy 29 (2004) 1825-1842 


The heat power Q\y includes the leakages through the side of the collector 
exposed to the sun and by the side in contact with the inside of the engine. The 

heat power discharged by the engine into the outside medium is 2 C oid- 

The fluxes of heat internal to the engine are represented with the help of the Stir¬ 
ling cycle imperfect exchange model, by Magot-Cuvru [9]. The heat power 
retained, Qhot, is divided into Q h, supplied to the driving medium, and 

(1 — a) | Greg |, supplied to the regenerator. The power discharged, Q co \d, includes the 
power Q c originating from the driving medium and (1 — a) | Q rQg | from the regener¬ 
ator. The regenerator efficiency a is such that (1 — a)|g reg | represents the heat orig¬ 
inating from the hot source, passing through the regenerator and discharged into 

the cold source, for each cycle. The driving medium, which retains gh and supplies 
• • 

Q c , also supplies the mechanical power indicated by W and the dissipated one 
\V f r i c . Simultaneously, during a cycle, it retains from and gives to the regenerator, 
the power g reg . 

The particular arrangements relative to the Stirling cycle must now be specified. 
During the intermediary evolutions (2-3) and (4-1), the pressure evolution is taken 
as being affine according to the volume, as it is linear for the Bouasse cycle. This 
hypothesis amounts to interpolating linearly, in the (p, V) diagram, the evolution 
between points 2 and 3, on the one hand, and points 4 and 1 on the other. It is 
simple to verify that this modelling respects condition (5) on regeneration, when 
the condition of the Stirling cycles — (S 2 — S\) = S 4 — S 3 has been imposed. This 
modelling was the one chosen to facilitate further calculations. 

Finally, dissipation is taken into account. The mechanical work dissipated 
through the forces of viscosity acting in the driving medium flow and the forces of 
solid friction acting on the mechanisms is taken into account by the efficiency rjmec, 
defined by 


bImec 


\§p&v 

- | 

W f rot | 


§p dV 



where \§p dV\ represents the mechanical work of the reversible cycle. 
3.2. Global efficiency of the solar heat engine 


The global efficiency of the solar heat engine rj g can be written as 

hg VcoW * hmoi' (?) 

The definition of solar collector efficiency (for example, [12: p. 411]) and Eq. (6) 
enable us to write rj C o\\ = 1 — 1 2ik,coii | / Qs • By explaining the heat power values with 
regard to the fluxes of heat = h {k • (T M - r amb ) and </> lk = h {k • (T abs - T amb ), it 
can be deduced that 



T-dmb 


7~amb 
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with 


T M — ^~amb + 


0s 

h\k 



for the absorber boundary temperature, obtained when (f) s = </> lk . 

To express the engine efficiency, it must be noted that 2hot is expressed, in the 
present case of the Stirling cycle, as 2hot = 2h + (1 — a)l2reg| ? and 2coid is given by 
Scold = Sc - (1 - a)lSreg| (see Fig. 3). Thus, Eq. (3) expressing J? mot is given as 


^7 mot 


2 h + 2c 


2h T (i 2reg 


With the hypothesis of isotherm evolution reversibility, 2h = T max • (S 4 — S 3 ) and 
2c = T m i n • (S 2 — S\) can be noted. First principles applied to the evolution (2-3) 

enables one to write 2reg = U 3 — U 2 + d V, which, with the hypothesis of the 
law of sure evolution, becomes 


Sreg = mCy{T max - T^ n ) + ■ ( V 3 - V 2 ). (10) 

In addition, the variations of entropy, in the case of the present driving medium, 
can be expressed as S 4 — S 3 = mrln(V/[/V 3 ) and S\ — S 2 = mr\n(V\/V 2 ) and this 
leads to 


^7 mot 


T Xn^-T InF 

1 max 111 t * 1 min 111 J/r 

V 3 * 2 


T, 


max 


Vi, 1 — a (1 — a) 


In- 1 - 


V 3 y — 1 


+ 


1 - 


V 2 

V 3 


-T, 


mm 


1 — a (1-a )fV 3 


7-1 


V 2 


-1 


( 11 ) 


Finally, taking the condition of Stirling cycles S 4 — S 3 = S\ — S 2 into account, 
by writing p = V 4 /V 3 = V\ / V 2 and p 0 = V 1 /V 3 , we obtain 

Tmaxkip T m [ n \np 

^7 mot 




max 


\ | 1 - « | (1 -ot)(p-p 0 ) 


y -1 


2 P 


- T, 


mm 


l-a (1 — «)(p — Po) 


y -1 


2 p 0 


( 12 ) 


The final stage consists in taking the irreversibility values into account. As they 
only appear due to dissipation and they only affect the work and not the exchan¬ 
ged heat values, it is sufficient to multiply by the mechanical efficiency rjmec- The 
efficiency can be finally expressed as 


^ 7 mot ^ 7 mec 


T max l up T m i n lnp 


T 

1 

r 

- 

1 

& 

- 

1 

sy 

1 

0 

T 

1- 

0 

1 

TT 

1 

1 

1_ 

1 max 

r^'y-V 2 p 

1 mm 

—1 

0 

<N 

1 

_1 


(13) 







































1834 


L.-C. Valdes / Renewable Energy 29 (2004) 1825-1842 


By substituting (8) and (13) in (7), we can express the global efficiency o g of the 
solar heat engine: 


0 * = n 


^abs — ^amb \ ^ ^ i \ 

1 7V E ) (Tmax hip T m i n mp) 

1 M d am b J 


mec 


p 


max 


\ | 1 ~ « | (i -V(p-po) 


y -1 


2 P 


- T, 


mm 


l-ct (1 — a)(p — p 0 ) 


7 - 1 


2 p 0 


(14) 


3.3. Heat power exchanged with the sources 


The value of engine power is entirely determined by the heat power values of 

Shot and Qcoid- The latter are imposed by the laws of heat exchange in the driving 
medium with the hot and cold sources. It is accepted that the heat exchanges are 
indeed represented by the laws of Newton, with coefficient heat transfer //hot and 
^coid? respectively, taken to be constant. 

The heat power supplied by the absorber is thus Q^ ot = A a b s • /Jhot(T a bs — T max ). 

By substituting this expression of ghot i n (6) and by explaining Q$ in an analogical 
way, we obtain 

(7m - T abs ) = ^(r abs - T max ). (15) 

h \k 

Likewise, the heat power lost by the engine is gcoid = -4 C oid * ^coid(T am b — T m i n ). 

In addition, relation (3) leads to gcoid = — (1 — ^ m ot)2hot- By using (15) to elimin¬ 
ate T max , we obtain 


(1 hmot) ' (Tm T a bs) 


^cold ' ^cold 
h\g ' ^abs 


(T m i n - T amb ) 



4. Optimised solar heat engines 

4.1. Principle of optimisation 

Optimising the dimensions of the solar heat engine requires notions of finite time 
thermodynamics, described by Feidt [12]. 

The problem being treated contains variables, values for which dimensioning 
occurs, in this case T max , T m i n , p, p 0 and a. The variables that are inde¬ 

pendent of dimensioning, /z hot , h co id, h^, 0 S , T amb and, according to (9), T M are 
considered as parameters. 

The problem consists in expressing, according to the parameters, the variables 
that maximise global efficiency, given by (14), under the constraints in (15) and 
(16). Studying the partial derivatives of rj g , established on the basis of (14), easily 
demonstrates that the efficiency is a monotonous function of the variables p, p 0 
and a. Optimising according to these values consists, therefore, in using the smal¬ 
lest value or the largest concretely permitted value, according to whether rj decrea¬ 
ses or increases with the considered variable. 
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The mathematical form of the constraint in (16) makes application of Lagrange’s 
method of multipliers difficult. For this reason, optimising according to T abs was 
studied directly which first required T max and T min to be eliminated from (14) using 
(15) and (16). The optimised dimensioning of the solar engine corresponds to the 
variables relative to the value of T abs which makes r\ g maximum. 

4.2. Expression of efficiency in conditions close to optimum 

For easier calculation, it is useful to introduce the variable 


X = 




and reason using the reduced temperatures of E max /EM and E^/Em- With these 
notations, (15), which can be written as 


F 


max 


T] 


M 


(1 +a 0 )X -a 0 , 


by writing 



0 o 


h\k 

^ch 


(19) 


supplies the relation giving E max /EM according to X. 

Relation (16) brings the factor 1 — rj mot into play. By replacing rjmot with the 
expression at (13), o g appears as a rational fraction in X having a degree-4 poly¬ 
nomial as the numerator as well as the denominator, which notably complicates 
the expression of the maximum. To get around this difficulty, we can suppose that 
the final result will not be significantly removed from the result obtained in [12: p 
333]; we can, therefore, admit at this stage in the calculation that 

1 — *7mot ~ v T a mb / E M . With this simplification, (16) can be written as 


F 


mm 


Em 


= —bob\X + (b 0 bi +b 2 ) 


( 20 ) 


by writing 


bo = 


h\k ' ^abs 
^cold ' ^cold 


^ ^1 1 ^7mot 

bi 


Lamb 


Ei 


( 21 ) 


Famb 


M 


Em 
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Using the notations in (21) and by writing 




< 


V 


K 

Co 

Cl 

C2 


n 


mec 


1 - b 2 

In p 

1 — ol (1 


7 - 1 
1 — a 

7 - 1 


+ 


(i 


«)(p - Po) 
2 P 

«KP ~ Po) 

2 po 


the expression of global efficiency rj g in (14) can also be written as 



Unax Unin 

rj =Kc 0 (l-X) - bM - Tm , (23) 

/ . \ J- max J- min 

(co + c\)— - c 2 -=— 

-L M 4 m 

and the elimination of the intermediary variables Unax /7m and Unin/Uvi using the 
relations in (18) and (20) finally supplies 


n 


g 


K( 1-X) 


(yiXjzlh) 

(a 2 X - p 2 ) ’ 


where 



fai = l+a 0 + bob\ 

\ = cio + bobi H- b 2 ’ 


and 





(1 + cio) H bob\ 
Co 


do + 


— bobi + 
Co 




were written. Using these notations, the engine efficiency rj mot can be written as 


^7mot dmec 


(a 2 X - p 2 ) ' 



The optimal temperature r abs can be found using Eq. (24). The extremum con¬ 
dition df] g /dX = 0 leads to the maximum in X, given by the relation 



P\ | P\ P 2 \ _ Q 

oq oq ot 2 J 



As Pi/oci > /? 2 /a 2 and 1 > j3 2 /a 2 , ^e root to be retained can be expressed as 





V i 

fi 2 , /?2 

1 5 

V ^2 \ 

/ oq 

a 2 a 2 
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and gives, transferred in (24), the optimal efficiency relative to this approximate 
calculation: 


*7g,opt 


= K 


0 C\ 

gl 2 


X 1- 



/ft 

ft 

V a 2 \ 

/ ai 

a 2 


1 - 


@2 Pi 


@2 | P 2 


a 2 V ai a 2 a 2 


Pi_ 

0 C\ 


1 - 


$2 Pi P 2 


a 2 V ai 


a 2 


( 30 ) 


By substituting in Eqs. (29), (18) and (20), we obtain 


Tabs = T] 


M 


1 _ /111 ^2 | P 2 


\ E max = Tm 


Emin = Em 


a 2 y oil 

(1 + Go) 


a 2 a 2 


L ^2 | fe 

a 2 Vai a 2 a 2 


a 0 


(31) 


-Mi 


i-k h-k + k) + {b(jb{ + b2) 

a 2 V ai a 2 a 2 


4.3. Calculation of performances 

Eqs. (13)—(16) of the model take into account the most significant imperfections: 
imperfect regenerator and heat transfers carried out with non-nil temperatures dif¬ 
ferences. Optimisation leading to the formulas in (30) and (31), based on the 
approximation 1 — rj mot « ^/T^/Tm, rests on the hypothesis that the engine 
reaches Carnot efficiency. Numerically, we can compensate for this hypothesis 
using an iterative calculation, based on a fixed point method. 

This iterative calculation assumes as an entry variable and 1 — rj mot as an exit 
variable. It is carried out for the value X = X opt given by (29). It consists in calcu¬ 
lating 1 — rjmou according to b\, using Eq. (27), and the variables oq, a 2 , j3 1 and f} 2 
depending on b\, having been previously calculated with the help of (25) and (26). 
The initial evaluation of b\ is carried out, in the context of the approximate optimi¬ 
sation, using Eq. (21). The numerical applications carried out show that the 
method converges very rapidly. 

The value X = X opt given by (29) contains the hypothesis 1 — rj mot « s/f im ^/T M 
and consequently, does not correspond to the maximum value for the global 
efficiency rj g . It is possible, without any risk of divergence, to apply the above iter¬ 
ative method for arbitrary values of X quite close to X opt and thus search, by suc¬ 
cessive approximations, for the value of X which makes the global efficiency rj g 
maximum. Such a calculation gives the exact optimal characteristics. 
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5. Example of a competitive solar heat engine 

5.1. Data 

The collecting area is fixed at S a bs = 1 m 2 . The solar heat engine functions with 
a solar flux of heat 0 S = 800 W/m 2 , in ambient atmosphere at the temperature 
Ta m b = 293 K. The choice of air as the driving medium leads us to adopt r = 
287 J/kg K and y = 1.4. 

By manufacturing the regenerator according to the model indicated in Section 

2 .2. effectiveness a close to 1 may be obtained by increasing the number of stages. 
Calculations, which will soon be published, demonstrate that the value of oc = 0.9 is 
a realistic one. The total rate of compression p 0 = V\ / V 3 cannot reach very high 
values in free-displacer engines. The value of po = 3 is compatible with the above 
value of a. The rate of compression p, linked to the duration of the intermediary 
phases, is fixed at 2.8. 

In alternative conventional engines, a significant proportion of the friction origi¬ 
nates from solid friction between the piston segments and the cylinder. In these 
engines, friction is unavoidable as it ensures air-tightness. It should be noted that 
power lost through friction increases with reaming and the power produced with 
capacity, so the piston-cylinder solution is only of value with sufficiently large 
engines of unitary power. The innovative solution in Section 2.2 completely elim¬ 
inates this inconvenience, in so much as friction is reduced to that created in the 
mechanical system of movement transformation and in the driving medium inter¬ 
nal flows. As all this friction diminishes with the rotational speed of the engine, it 
is possible to fix it at a low value, giving very good mechanical efficiency. The 
value of Omec = 0.9, with the preceding values and global efficiency of approxi¬ 
mately 0.15, leads to dissipated power of the order of 12 W, which is realistic, even 
if slightly overestimated. 

The driving medium thicknesses on the absorber and cooler are of the order of a 
few millimeters, as a result of the surfaces of largest heat exchange in view of the 
volumes, of the order of liters, generated by the free piston. The thermal resistance 
of the absorber being essentially due to the thickness of the air, we are led to adopt 
Miot = 50W/m 2 K using a simple calculation of transference by conduction. 
Regarding the cooler, the resistance due to the natural convection on the fan is 
added to this thermal resistance. The value of /z co id = 10 W/m 2 K is adopted. The 
load piston, located on the fan side, reduces the exchange area ^4 co id- We accept 
A coU = 0.85 m 2 . 

5.2. Results and comments 

The methods of calculation in Section 4.3, giving the optimal values of the 
absorber temperature r abs , global efficiency rj g and the temperatures T max and 
T m in, are applied to the numerical data above. We chose to study the variation of 
these variables according to /z ]k , as it is the value that can be directly interpreted in 
technical terms, and therefore, economic competitiveness. 
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Table 1 


Evolution of values with h\ k . Approximate calculation 


Aik [W/m 2 °C] 

5.46 

2.73 

1.37 

0.91 

0.68 

0.55 

0.46 

^s/^lkTamb 

0.5 

1 

2 

3 

4 

5 

6 

T m [°C] 

166.5 

313.0 

603.9 

899.1 

1196.5 

1474.5 

1759.1 

7/bs [°C] 

117.9 

180.4 

273.7 

349.4 

415.1 

470.1 

521.6 

r max [°C] 

112.6 

173.1 

264.7 

339.5 

404.5 

459.0 

510.3 

T r ni„ [°C] 

45.5 

50.1 

50.8 

49.4 

47.9 

46.6 

45.4 

*7g 

0.050 

0.106 

0.186 

0.240 

0.280 

0.309 

0.334 

h /mot 

0.151 

0.234 

0.328 

0.384 

0.422 

0.448 

0.469 

^Carnot 

0.250 

0.354 

0.465 

0.529 

0.574 

0.606 

0.631 


Table 1 gives the numerical results obtained by applying Eqs. (30) and (31), rela¬ 
tive to approximate optimisation. Table 2 gives the results obtained using the iter¬ 
ative method applied with a calculated absorber temperature of T a b s , according to 
the approximate optimisation, using Eq. (29). Table 3 demonstrates several exact 
calculations obtained using the method of optimisation through successive approx¬ 
imations. Finally, Fig. 4 represents the evolution of working efficiency and the loss 
coefficient of h\y, and Fig. 5, that of temperatures. 

Figs. 4 and 5 comply with the fact that the absorber temperature and efficiency 
increase as the heat leakages diminish. For the values of h\y > 0.46 W/m 2 K, the 
obtained values of T ahs correspond to the electromagnetic radiation of the absorber 
in the visible domain. As the real heat leakages are, therefore, no longer compat¬ 
ible with those anticipated, taken to be without loss by radiation, the evolution is 
only represented up to h\y > 0.46 W/m 2 K. 

The numerical results reveal that the efficiency of the Stirling engine studied is 
always between 60% and 74% of the Carnot efficiency, whereas the global efficiency 
is, due to the heat leakages from the solar collector, between 20% and 53% of the 
Carnot efficiency. 

The values of which enable solar cell efficiencies of 0.12, 0.17 and 0.24 to be 
reached, are 2.5, 1.55 and 0.9 W/m K, respectively. The efficiencies of 0.12 and 
0.17 for ordinary solar cells can, therefore, be attained by solar engines of the 
type studied, when they are equipped with relatively usual solar collectors, whereas 
the efficiency of 0.24 requires flat plate collectors giving an improved, but not 


Table 2 


Evolution of values with h\ k . Iterative calculation to the approximate optimum 


A, k [W/m 2 °C] 

5.46 

2.73 

1.37 

0.91 

0.68 

0.55 

0.46 

T m [°C] 

166.5 

313.0 

603.9 

899.1 

1196.5 

1474.5 

1759.1 

T abs [°C] 

117.9 

180.4 

273.7 

349.4 

415.1 

470.1 

521.6 

T m , x [°C] 

112.6 

173.1 

264.7 

339.5 

404.5 

459.0 

510.3 

T min [°C] 

46.6 

52.8 

56.2 

56.8 

56.7 

56.5 

56.2 

*7g 

0.049 

0.104 

0.181 

0.235 

0.274 

0.303 

0.327 

7/mot 

0.149 

0.229 

0.321 

0.375 

0.413 

0.439 

0.459 

Carnot 

0.250 

0.354 

0.465 

0.529 

0.574 

0.606 

0.631 
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Table 3 

Evolution of values with h {k . Precise calculation 


h ik [W/m 2 °C] 

5.46 

0.91 

0.46 

T m [°C] 

166.5 

899.1 

1759.1 

T „ bs [°C] 

120.2 

367.6 

553.1 

r max [°C] 

115.1 

357.9 

541.9 

T t ni„ [°C] 

45.1 

54.8 

54.5 


0.049 

0.235 

0.327 

^1 mot 

0.156 

0.389 

0.472 

Carnot 

0.255 

0.543 

0.645 


exceptional, performance. Finally, if we consider, as the numerical example given 
by Pasquetti [13: p. 7] leads us to believe, that the value h\y = 0.50 W/m 2 K is 
technically acceptable, then the graph in Fig. 5 indicates that a global efficiency of 
0.31 may be reached. This greatly exceeds the maximal efficiency of 0.24 for solar 
cells. 

To finish, we can verify that the approximate calculation supplies efficiency 
values less than 4% of the exact value and the iterative calculation, less than 0.01%. 
Approximate optimisation leads us to underestimate the temperature: from 3 to 30 °C 
for the absorber temperatures and maximal for the driving medium, from 1 and 
10 °C for the driving medium minimum temperatures. 



Fig. 4. Evolution of efficiency with heat leakages. Example of Section 5.1. 
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Fig. 5. Evolution of temperatures with heat leakages. Example of Section 5.1. 


5.3. Evaluation ofkWh costing 

The production of a first prototype with a collector area of 1 nr at the GREEn 
has enabled the cost of a solar heat engine according to the model in Section 2.2 to 
be estimated. On the basis of commercial elements, the cost of the solar engine in 
industrial production was evaluated at €525. With a lifespan of 20 years, with an 

r\ 

annual average of 1760 kW h/nr year, and with 0.15 efficiency, including the 
efficiency value of 0.5 for energy storage and destocking in an accumulator battery, 
we obtain a cost price of €0.099/kW h, perfectly comparable with that of the kW h 
produced by power stations. 

6. Conclusion 

Replacing the piston-cylinder system with metallic bellows solves the problem of 
air-tightness without energy dissipation. This innovation opens the possibility of 
transforming heat into work, with good real efficiency, using alternative heat 
engines of low power. 

The solar heat engine studied, of low unitary mechanical power (100-3000 W), 
integrates this innovation. To guarantee its economic competitiveness, it includes a 
flat plate solar collector and a free-displacer Stirling engine. Its competitiveness is 
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not, however, guaranteed unless its dimensioning is optimised. To this end, notions 
of finite time thermodynamics are applied to a mathematical model, which takes 
into account the most significant penalising factors: differences in dynamic tem¬ 
peratures between the sources and the driving medium, rate of compression and 
regenerator imperfections. 

This study demonstrates that the above innovation enables the 0.24 efficiency 
value of the best solar cells to be exceeded. Strictly from the point of view of 
efficiency, such systems are serious competitors to solar cells. In addition, the tech¬ 
nological production of these machines enables a moderate cost price. Approxi¬ 
mate figuring, which gives a kW h cost in the order of €0.01, allows us to think 
that this system is equally competitive with regard to power stations. 

These elements should enable the heat option of exploiting solar electromagnetic 
radiation to become a possibility. 
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